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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2646 

INVISCID  FLOW  ABOUT  AIRFOILS  AT  HIGH  SUPERSONIC  SPEEDS 
By  A.  J.  Eggers,  Jr*,  and  Clarence  A.  Syvertson 

SUMMARY 


Steady  flow  about  curved  airfoils  is  investigated  analytically  at 
high  supersonic  speeds.  Assuming  air  "behaves  as  an  ideal  diatomic  gas, 
it  is  found  that  small  pressure  disturbances  emanating  from  the  surface 
of  an  airfoil  are  almost  completely  absorbed  in  the  leading-edge  shock 
wave  (or  a  shock  wave  emanating  from  any  other  location  on  the  surface), 
provided  the  flow  deflection  angles  are  not  too  close  to  those  corre¬ 
sponding  to  shock  detachment.  This  result  is  found  to  be  essentially 
independent  of  Mach  number.  As  a  consequence,  it  is  shown  that  within 
the  limitations  of  the  assumption  of  ideal  gas  flow,  the  shock -expans ion 
method  may  be  used  with  good  accuracy  to  predict  pressure  distributions 
on  curved  airfoils  at  arbitrarily  high  Mach  numbers.  This  observation 
is  verified  with  the  aid  of  the  method  of  characteristics  applied  to  a 
10-percent-thick .biconvex  airfoil  at  0°  angle  of  attack.  It  is  further 
shown  that  the  shock-expansion  method  can  be  easily  employed  to  con¬ 
struct  the  entire  flow  field  about  a  curved  airfoil,  accounting  for 
shock-wave  curvature  and  resulting  entropy  gradients  in  the  flow. 

An  approximation  to  the  shock-expansion  method  for  thin  airfoils 
at  high  Mach  numbers  is  also  investigated,  and  is  found  to  yield  pres¬ 
sure  distributions  in  error  by  less  than  10  percent  at  Mach  numbers 
above  3  and  flow  deflection  angles  up  to  25°.  This  slender -airfoil 
method  is  relatively  simple  in  form  and  thus  may  prove  useful  for  some 
engineering  purposes.  To  this  end,  tables  are  presented  to  facilitate 
its  use. 

Effects  of  caloric  imperfections  of  air  manifest  in  disturbed  flow 
fields  at  high  Mach  numbers  are  investigated,  particular  attention  being 
given  to  the  reduction  of  the  ratio  of  specific  heats  from  1.4  toward 
1.0.  So  long  as  this  ratio  does  not  decrease  appreciably  below  1.3*  it 
is  indicated  that  the  shock-expansion  method,  generalized  to  include 
effects  of  these  imperfections,  should  be  substantially  as  accurate  as 
for  ideal  gas  flows.  This  point  is  checked  by  comparing  pressure  distri¬ 
butions  predicted  by  the  generalized  shock-expansion  method  and  a  gen¬ 
eralized  method  of  characteristics.  Both  methods  are  employed  in  forms 
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applicable  for  local  air  temperatures  up  to  about  5000°  Rahkine,  corre¬ 
sponding,  for  slender  airfoils,  to  Macb  numbers  up  to  the  order  of  the 
so-called  escape  Mach  number.  Caloric  imperfections  caused  reductions 
in  the  pressure  coefficients  below  those  predicted  for  flows  of  an 
ideal  gas.  In  turn,  there  is  a  general  reduction  in  force  and  moment 
coefficients  up  to  10  percent. 

The  slender-airfoil  method  is  modified  to  employ  an  average  value 
of  the  ratio  of  specific  heats  for  a  particular  flow  field.  This  sim- 
plified  method  has  essentially  the  same  accuracy  for  imperfect  gas 
flows  as  its  counterpart  has  for  ideal  gas  flows. 

An  approximate  flow  analysis  is  made  at  extremely  high  Mach  numbers 
where  it  is  indicated  that  the  ratio  of  specific  heats  may  approach 
close  to  1.  It  is  found  that  the  shock -expans ion  method  may  be  in  con-, 
siderable  error  as  disturbances  incident  on  the  leading-edge  shock  wave 
are  no  longer  largely  absorbed  in  the  wave.  In  this  case,  however,  the 
Busemann  method  for  the  limit  of  infinite  free-stream  Mach  number  and 
specific  heat  ratio  of  1  appears  to  apply  with  reasonable  accuracy. 


INTRODUCTION 


Small  disturbance,  potential— flow  theories  have  been  employed 
widely,  and  for  the  most  part  successfully,  for  predicting  the  pressures 
(and  velocities)  at  the  surface  of  an  airfoil  in  steady  motion  at  low 
supersonic  speeds.  Thus  the  linear  theory  of  Ackeret  (reference  l)  has 
proven  particularly  useful  in  studying  the  flow  about  relatively  thin, 
sharp-nosed  airfoils  at  small  angles  of  attack,  while  the  second-order 
theory  of  Busemann  (reference  2)  has  found  application  when  thicker 
airfoils  at  larger  angles  of  attack  were  under  consideration.  At  high 
free-stream  Mach  numbers  the  range  of  applicability  of  any  potential 
theory  is  seriously  limited,  however,  due  to  the  production  of  strong 
shocks  by  even  the  relatively  small  flow  deflections  caused  by  thin 
airfoils.  The  assumption  of  potential  flow  is  invalidated,  of  course, 
by  the  pronounced  entropy  rises  occurring  through  these  shocks. 

This  limitation  on  potential  theories  was  early  recognized  and  led 
to  the  adoption  (see  reference  3)  of  what  is  now  commonly  called  the 
shock— expansion  method.  The  latter  method  derives  its  advantage  over 
potential  theories  principally  by  accounting  for  the  entropy  rise 
through  the  oblique  shock  emanating  from  the  leading  edge  of  a  sharp¬ 
nosed  airfoil.  Consequently,  so  long  as  the  disturbed  air  behaves 
essentially  like  an  ideal  gas,  and  so  long  as  entropy  gradients  normal 
to  the  streamlines  (due  to  curvature  of  the  surface)  do  not  signifi¬ 
cantly  influence  flow  at  the  surface,  the  shock-expansion  theory  should 
predict  the  pressures  at  the  surface  of  an  airfoil  with  good  accuracy  - 
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it  is  tacitly  assumed,  of  course,  that  the  flow  velocity  is  everywhere 
supersonic,  and  that  the  Reynolds  number  of  the  flow  is  sufficiently 
large  to  minimize  viscous  effects  on  surface  pressures. 

The  departure  of  the  behavior  of  air  from  that  of  an  ideal  gas  at 
the  temperatures  encountered  in  flight  at  high  supersonic  speeds  has 
been  the  subject  of  some  investigation  in  the  case  of  flows  through 
oblique  shock  waves.  In  reference  4,  the  effects  of  thermal  and  caloric 
imperfections  on  the  pressure  rise  across  an  oblique  shock  wave  was 
investigated  at  sea-level  Mach  numbers  of  10  and  20  and  it  was  found 
that  these  effects  decreased  the  rise  by  less  than  5  percent  for  maximum 
temperatures  up  to  3000°  R  (corresponding  to  flow  deflection  angles  up 
to  24°) .  This  decrease  was  found  to  be  due  almost  entirely  to  caloric 
imperfections  or  changes  in  vibrational  heat  capacities  of  the  air  pass¬ 
ing  through  the  shock  wave.  The  changes  in  temperature  and  density  of 
the  air  passing  through  the  wave  were  affected  to  a  considerably  greater 
extent.  Subsequently,  an  investigation  was  carried  out  by  Ivey  and 
Cline  up  to  Mach  numbers  as  high  as  100  (reference  5)  using  the  results 
for  normal  shock  waves  obtained  by  Bethe  and  Teller  considering  effects 
of  dissociation  (reference  6) .  As  would  be  expected,  the  pressures 
were  found  to  be  affected  to  a  somewhat  greater  extent  at  the  higher 
Mach  numbers. 

The  extent  to  which  flow  in  the  region  of  the  leading  edge  of  an 
airfoil  departs  from  the  simple  Prandtl -Meyer  type  has  also  been  inves¬ 
tigated  at  high  supersonic  airspeeds.  If  the  surface  is  curved,  for 
example,  to  give  an  expanding  flow  downstream  of  the  leading  edge, 
expansion  waves  from  the  surface  will  interact  with  the  nose  shock  wave, 
thereby  curving  it  and  yielding  a  nonisentropic  flow  field.  This  flow 
field  may  be  characterized  not  only  by  disturbances  emanating  from  the 
surface  but  also  by  disturbances  reflecting  to  some  extent  from  the 
shock  wave  back  toward  the  surface.  The  manner  in  which  these  phenomena 
dictate  shock-wave  curvature  and  surface  pressure  gradient  in  ideal  gas 
flows  at  the  leading  edge  has  been  treated  by  Crocco  (reference  7)  and 
more  recently  by  Schaefer  (reference  8),  Munk  and  Prim  (reference  9)> 
and  others.  In  the  cases  considered  by  Munk  and  Prim  it  was  found  that 
surface  pressure  gradients  were  less  (in  absolute  value)  than  those 
obtained  assuming  Prandtl -Meyer  flow  at  the  higher  Mach  numbers  (i.e., 
Mach  numbers  greater  than  about  3)  although  generally  by  no  more  than 
about  10  percent.  Since  curved  airfoils  are  likely  to  be  of  fundamen¬ 
tal  interest  at  high  flight  speeds  (see,  e.g.,  reference  10) ,  these 
phenomena  would  appear  to  merit  further  investigation,  particularly  as 
regards  their  influence  on  the  whole  flow  field.  In  addition  it  would 
appear  desirable  to  consider  effects  of  gaseous  imperfections  through¬ 
out  the  field. 
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Such  an  investigation  has  therefore  been  undertaken  in  the  present 
report  uging  the  method  of  characteristics  to  accurately  obtain  flow 
fields,  and  as  a  basis  for  obtaining  the  more  approximate  methods  of 
analysis.  The  method  is  employed  in  a  generalized  form  which  allows 
caloric  imperfections  as  well  as  entropy  gradients  in  the  flow  to  be 
considered  at  tenperatures  up  to  the  order  of  5000°  R  -  thermal  imper¬ 
fections  are  neglected  as  being  unimportant  in  atmospheric  air  flows 
(see  reference  4)*  A  10-percent-thick  biconvex  airfoil  is  treated  at 
Mach  numbers  from  3*5  to  infinity,  and  the  results  are  compared  with 
the  predictions  of  the  shock -expans ion  method,  including  a  simplified 
form  of  the  method  applicable  to  slender  airfoils  at  high  Mach  numbers, 
and  a  generalized  form  of  the  method  including  effects  of  caloric 
imperfections. 


SYMBOLS 


a 


c 


'l* 


Cd 


c 


p 


cv 


M 

P 

q. 

R 


local  speed  of  sound,  feet  per  second 
chord,  feet 

characteristic  coordinates  (Cx  positively  inclined  and  C2 
negatively  inclined  with  respect  to  the  local  velocity 
vector) 

section  drag  coefficient 
section  lift  coefficient 

section  moment  coefficient  (moment  taken  about  leading  edge) 


pressure  coefficient 


Po 


<lo 


specific  heat  at  constant  pressure,  foot-pounds  per  slug  °R 
specific  heat  at  constant  volume,  foot-pounds  per  slug  °R 
Mach  number  (ratio  of  local  velocity  to  local  speed  of  sound) 
static  pressure,  pounds  per  square  foot 
dynamic  pressure,  pounds  per  square  foot 
gas  constant,  foot-pounds  per  slug  °R 
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s,n 

T 

t 

V 

x,y 

a 

P 

7 

6 

6 

P 

CT 

U) 

o 

A,B,  1 
C  ,D. • • J 

i 

N 

S 


rectangular  coordinates  (in  streamline  direction  and  normal 
to  streamline  direction,  respectively) 

temperature,  °R 
time,  seconds 

resultant  velocity,  feet  per  second 
rectangular  coordinates 

angle  of  attack,  radians  unless  otherwise  specified 
Mach  angle,  arc  sine  (h)  ,  radians 


ratio  of  specific  heats  l  — 

\cv 

(Average  value  of  y  is  7a-) 

flow  deflection  angle,  radians  unless  otherwise  specified 

molecular  vibrational  energy  constant,  °R  (5500°  R  for  air) 

mass  density,  slugs  per  cubic  foot 

shock -wave  angle,  radians 

ray  angle  for  Prandtl-Meyer  flow,  radians 

Subscripts 


free -stream  conditions 

conditions  at  different  points  in  flow  field 
ideal  gas  quantities 

conditions  just  downstream  of  shock  wave 
conditions  on  streamline 

Superscript 


vector  quantities 
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DEVELOPMENT  OF  METHODS  OF  ANALYSIS 
Method  of  Characteristics 


Two-dimensional  rotational  supersonic  flows  have  been  treated  by 
numerous  authors  with  the  aid  of  the  method  of  characteristics,  and 
various  adaptations  of  the  method  have  been  found  which  are  especially 
suited  for  studying  particular  types  of  such  flows.  In  the  case  of 
steady  flows  in  which  atmospheric  air  does  not  behave  as  an  ideal  dia¬ 
tomic  gas,  a  very  familiar  and  simple  form  of  the  compatibility  equa¬ 
tions  may  be  employed.  To  illustrate,  consider  the  Euler  equation 


P 


dV 

dt 


-  grad  p 


(1) 


the  continuity  equation 


div  (pV)  =0  (2) 

and  the  equation  for  the  speed  of  sound  (evaluated  at  constant  entropy) 


Rewriting  equations  (l)  and  (2)  in  the  form  of  partial  differential 
equations  and  transforming  the  resulting  expressions  to  the  character¬ 
istic  or  Ci,  C2  coordinate  system,  there  is  obtained,  upon  combination 
with  equation  (3),  the  following  relations  for  steady  flow: 


and 


ctg  P  /  dp  dp  \  f  &  ,  &  \ 

pV2  V^Ci  dc2/  \5Ci  bC2J 


ctg  3  /  dp  [  dp  \  &  d8_\ 

PV2  VdCi  dc2/  \  dCi  dc £J 


=  0 


=  0 


W 


(5) 


A  simple  addition  or  subtraction  of  equations  (4)  and  (5)  then  yields 
the  compatibility  equations  (see,  e.g.,  reference  11) 


dp 

dCi 


=  —  p  V2  tan 


(6) 
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and 

=  pV2  tan  p  (7) 

8c2  6c2 

Now,  in  reference  4  both  caloric  and  thermal  imperfections  of  air  were 
considered  and  it  was  found  that  the  latter  imperfections1  have  a  neg¬ 
ligible  effect  on  shock  processes  in  atmospheric  air.  It  may  easily  be 
shown  that  this  conclusion  also  applies  to  expansion  processes ,  and  for 
this  reason  caloric  imperfections,  only,  are  considered  in  detail  in 
the  present  paper.  These  imperfections  become  significant  in  air  at 
temperatures  greater  than  about  800°  R  and  first  manifest  themselves  as 
changes  in  the  vibrational  heat  capacities  with  temperature.  Thus,  the 
specific  heats,  Cp  and  cv,  and  their  ratio,  7,  for  the  gas  also  change. 
The  equation  of  state  remains,  however, 

P  =  pRT  (8) 

\ 

and  the  specific  heats  are  still  related  to  the  gas  constant  by  the 
'expression 

Cp  -  cv  =  R  (9) 


Furthermore,  it  readily  follows  from  the  differential  energy  equation 
and  these  expressions  that  the  speed  of  sound  is  given  by  the  simple 
relation 

a2  =  7RT  (10) 

Combining  equations  (8)  and  (10)  and  noting  that  sin  3  =  a/v  there  is 
then  ohta-ined 


TP 

sin2  p 


(11) 


Hence,  on  combining  this  equation  with  equations  (6)  and  (7),  it  is 
apparent  that  the  familiar  compatibility  equations 


8p  _  -2tp  8s 
8cx  sin  2  3  dCj 


(12) 


4 


Thermal  imperfections  usually  appear  in  the  form  of  intermolecular 
forces  and  molecular -size  effects,  and  may  be  accounted  for  with 
additional  terms  in  the  equation  of  state. 
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and 


5p _ 2  TP  5s 

Sc2  sin  2  p  Sc2 


(13) 


also  hold  for  the  more  general  type  of  flow  under  consideration.  These 
equations  are  basic,  of  course,  to  two-dimensional  characteristics 
theory,  and,  as  will  he  shown  later,  form  a  convenient  starting  point 
for  developing  simpler  theories  of  two-dimensional  supersonic  flow. 


In  order  to  apply  equations  (12)  and  (13),  it  is  evident  that  the 
manner  in  which  y  and  p  or  M  are  connected  to  p  or  8  must  he  known. 
Relations  implicitly  connecting  these  variables  at  temperatures  up  to 
the  order  of  5000°  R  may  he  readily  obtained  from  the  results  of  refer¬ 
ence  4  by  simply  eliminating  the  terms  therein  accounting  for  thermal 
imperfections.  Thus  we  have  as  a  function  of  the  local  static  tempera¬ 
ture  and  free-stream  conditions 


For  isentropic  flow  along  a  streamline,  the  pressure  is  related  to  the 
temperature  by  the  expression 


where 


(16) 


+  (T)  = 


(17) 
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If  there  is  a  shock  wave  in  the  flow,2  in  particular  a  nose  or  leading- 
edge  shock,  then  the  following  additional  relations  obtained  with 
equations  (8),  (10),  and  (15)  and  the  conditions  for  continuity  of  flow 
and  conservation  of  momentum  along  a  streamline  through  the  shock  are 
also  required: 


% 


T, 


(l+7#%2)""  ^  (l+7oMo2) 


J 


(l+7NMN2)-  —  (1+roMo2) 


■N 


k  — 

% 

. 

(18) 


and 


sin2o 


7nTnmn2  ^  i 

r0T0M02 


tan  6  =  - 

tan  a 


1 

7oMoa 

(PN/p0)-l 


(19) 


(20) 


Using  the  local  static  temperature  as  a  parameter,  the  term 
27p/sin  2  |3  in  equations  (12)  and  (13)  may  now  he  evaluated  with  equa¬ 
tions  (l4)  through  (17).  Equations  (l8)  through  (20)  define  the  initial 
conditions  downstream  of  a  leading-edge  or  other  shock  wave  in  the  flow 
field.  Thus,  equations  (12)  through  (20)  provide  all  the  information 
necessary  to  calculate  the  flow  about  an  airfoil  by  means  of  the  method 
of  characteristics.  As  described  in  detail  in  appendix  A,  the  calcula¬ 
tion  is  of  three  general  types;  namely,  (l)  calculation  of  conditions  at 
a  point  in  the  flow  field  between  the  shock  and  the  surface;  (2)  calcula 
tion  of  conditions  at  a  point  on  the  surface;  and  (3)  calculation  of 
conditions  at  a  point  just  downstream  of  the  shock.  Case  (l)  entails 
the  use  of  both  compatibility  equations,  while  case  (2)  entails  the 
use  of  the  compatibility  equation  for  a  second-family  characteristic 


2If  there  are  no  shock  waves,  then  the  subscript  N  in  equation  (l6) 
can,  of  course,  be  replaced  with  the  subscript  o.  _ 
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line  in  combination  with  the  equation  of  the  airfoil  surface,  and 
case  (3)  involves  the  compatibility  equation  for  a  first-family  line 
in  combination  with  the  oblique  shock  equations.  With  the  aid  of  the  * 

three  general  types  of  calculations  the  entire  flow  field  about  an 
airfoil  can  be  built  up  numerically  using  a  computing  procedure  working 
from  the  leading  edge  downstream.  In  cases  where  changes  in  the  vibra¬ 
tional  heat  capacities  with  temperature  are  neglected,  the  calculations 
are  of  course  simplified  since  y  of  the  gas  can  be  considered  con¬ 
stant,  and  temperature ,  pressure,  and  density  ratios  are  simply  the 
ideal  gas  functions  of  Mach  number. 


Shock-Expansion  Method 


General.-  This  method  of  calculating  supersonic  flow  of  an  ideal 
gas  at  the  surface  of  an  airfoil  is  well  known,  entailing  simply  the 
calculation  of  flow  at  the  nose  with  the  oblique  shock  equations  and 
flow  downstream  of  the  nose  with  the  Prandtl-Meyer  equations.  Deter¬ 
mination  of  airfoil  characteristics  in  this  manner  requires  a  small 
amount  of  time,  of  course,  compared  to  that  involved  when  the  method  of 
characteristics  is  used,  hence  the  advantage  of  the  former  method.  The 
questions  arise,  however,  as  to  exactly  what  the  simplifying  assumptions 
underlying  the  shock-expansion  method  are,  and  what  form  the  method 
takes  (for  calculative  purposes)  when  the  gas  displays  varying  vibra¬ 
tional  heat  capacities. 

* 

The  matter  of  simplifying  assumptions  may  perhaps  best  be  con¬ 
sidered  by  employing  equations  (12)  and  (13),  the  basic  compatibility  4 

equations.  If  these  expressions  are  resolved  into  the  streamline 
direction  and  combined,  noting  that 


=  - 1 _  ($L+ 

ds  2  cos  0  VScx  dc2/ 

and 

d§>_  =  1  f  _d6_  +  ds  \ 

cte  2  cos  0  \  dCj  dC2/ 
there  is  then  obtained  the  relation 


(21) 


(22) 


4 


Sp 

ds 


sin  2  p  ds 


(23) 
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defining  the  gradient  of  p  along  s.  If  flow  along  streamlines  down¬ 
stream  of  the  nose  is  of  the  simple  Prandtl -Meyer  type,  however,  we 
have 


dp  _  2yp  dS 
ds  sin  2  p  6s 


(24) 


Hence  it  is  evident  that  the  requirement  for  this  type  of  flow  is 


cfc/dCi 

bsfoCz 


<  <  1 


(25) 


Equation  (25)  is,  of  course,  simply  an  approximate  statement  of  a  well- 
known  property  of  Prandtl -Meyer  flows;  namely,  that  flow  inclination 
angles  are  essentially  constant  along  first-family  Mach  lines.  It  fol¬ 
lows  from  equation  (12)  that  if  equation  (25)  holds,  then  the  pressures 
will  also  he  essentially  constant  along  these  lines.  It  does  not  fol¬ 
low,  however,  that  the  Mach  number  will  he  constant,  or  for  that  matter 
that  the  first-family  characteristic  lines  will  he  straight  (as  is  the 
case  for  isentropic  expansion  flows  about  a  corner).  In  fact,  it  may 
easily  he  shown  that  the  Mach  number  gradient  along  Ci  is  proportional 
to  the  local  entropy  gradient  normal  to  the  streamlines,  and  that  the 
Ci  lines  are  curved  according  to  the  change  in  M.  Thus  we  see  that 
there  is  really  only  one  basic  assumption  underlying  the  shock -expansion 
method;  namely,  disturbances  incident  on  the  nose  shock  (or  for  that 
matter  any  other  shock)  are  consumed  almost  entirely  in  changing  the 
direction  of  the  shock.3  Within  the  limitations  of  this  assumption  it 
is  evident  that  the  method  provides  a  relatively  simple  means  for  cal¬ 
culating  the  whole  flow  field  about  an  airfoil,  including  effects  of 
shock -wave  curvature  (see  appendix  B) .  In  general,  of  course,  the 
validity  of  this  assumption  can  only  be  checked  by  comparison  of  calcu¬ 
lations  using  this  method  with  those  using  the  method  of  characteristics. 

The  shock -expans ion  method  for  a  calorically  imperfect  diatomic  gas 
is  readily  deduced  from  the  equations  previously  obtained.  For  example, 
flow  conditions  at  the  leading  edge  of  an  airfoil  can  be  evaluated  with 
the  oblique  shock -wave  expressions  (equations  (l8)  through  (20))  and  the 
expression  for  conservation  of  energy  (equation  (15))*  The  variation 
of  flow  inclination  angle  with  pressure  along  the  surface  is  then 


3It  is  interesting  to  note  that  the  assumption  of  Thomas  (reference  12) 
that  pressure  is  a  function  only  of  flow  deflection  angle  and  entropy 
is  equivalent  to  this  assumption.  It  follows,  of  course,  that  the 
most  general  solution  obtainable  with  Thomas's  series  representation 
of  the  pressure  is  that  given  by  the  shock -expans ion  method. 
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obtained  by  graphically  integrating  equation  (24) ;  namely, 

r 

8s  -  8B  =  ^P»  dp  (26) 

ps 

where  the  variables  7,  p,  and  3  are  evaluated  using  equations  (l4) 
through  (17),  employing  the  static  temperature  as  a  parameter.  When 
extreme  accuracy  is  not  essential,  this  rather  tedious  calculation  can 
he  avoided  and  a  relatively  simple  algebraic  solution  of  the  flow  down¬ 
stream  of  the  nose  can  be  employed.4  The  details  of  this  solution  are 
presented  in  appendix  C.  In  the  special  case  of  flow  at  high  supersonic 
speeds  about  slender  airfoils,  the  whole  calculation  becomes  particularly 
simple  and  warrants  special  attention. 

Slender  airfoils » -  If  it  is  assumed  that  the  local  surface  slopes 
are  small  compared  to  1  and  in  addition  that  the  free -stream  Mach  number 
is  large  compared  to  1,  it  follows  that  0  and  3  are  everywhere  small 
compared  to  1.  In  this  case  equation  (24)  takes  on  the  approximate  form 


d5 


7PM 


(27) 


Furthermore,  if  it  is  assumed  that  7  is  constant  at  an  average  value 
7a  for  a  particular  flow  field  (this  assumption  appears  reasonable 
since  in  the  temperature  range  up  to  5000°  R  the  change  in  7  is  less 
than  10  percent  as  shown  in  reference  4) ,  then  the  Mach  number  and 
pressure  may  be  related  by  the  simple  expression 


Equations  (27)  and  (28)  combine  to  yield  the  differential  equation 


4 


1 

&  % 


=  as 


(29) 


which  readily  integrates  (between  N  and  S)  to  the  form 

4  The  tabulated  results  of  Noyes  (reference  13)  may  also  prove  useful  in 
this  case  for  Mach  numbers  up  to  3. _ 
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(30) 


(31) 


(32) 


there  is  obtained  from  the  oblique  shock  equations,  simplified  to  con 
form  with  this  analysis. 


f  (Mq8n)  = 


Mq2<Tn2  -  1 


and 


g(Mo8N)  = 


gyaMo2^2  ~  (/a"1) 
7*  +  1 


vhere 


Mq^N 


=  — i  M05N  +  J 1  +^1^1  MqSn) 


(33) 


(3*0 


(35) 


With  equations  (30)  through  (35)  the  pressures  on  the  surface  of  an 
airfoil  may  easily  he  obtained.  In  terms  of  pressure  coefficient  we 
have 


(36) 


or 


|  g(M05N) 


l-f(Mo5N)(h  -  ^ 


fls_ 

ya-1  -1 


(37) 
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The  advantage  of  these  slender-airfoil  expressions  lies,  of  course,  in 
their  relative  simplicity  and  thus  the  ease  of  calculation  which  is 
inherent  to  them.  It  may  he  noted  in  this  regard  that  the  functions 
f ( Mq8 jj- )  and  g can  he  calculated  once  and  for  all  with  equations 
(33),  (3*0,  and  (35),  provided  the  variation  of  ya  with  is  known. 

This  calculation  has  heen  carried  out  for  a  constant  value  of  y  equal 
to  1.4,  and  average  values  of  y  assuming  T0  =  500°  R.5  The  results 
are  presented  in  table  I. 

It  should  also  he  noted  that  the  slender-airfoil  expressions  of 
the  shock-expansion  method  satisfy  the  hypersonic  similarity  law  for 
airfoils  first  deduced  by  Tsien  (reference  l4) A  necessary  condition 
for  the  validity  of  these  expressions  is  thus  satisfied;  however,  the 
accuracy  of  the  shock-expansion  method,  whether  for  slender  airfoils  or 
otherwise,  remains  to  he  investigated.  Such  an  investigation  is  now 
undertaken  with  the  aid  of  the  method  of  characteristics. 


INVESTIGATION  OF  FLOW  ABOUT  AIRFOIIS 
AND  DISCUSSION  OF  RESULTS 


This  study  is  divided  into  two  parts:  first,  a  consideration  of  the 
effects  of  Mach  number  assuming  air  behaves  as  an  ideal  diatomic  gas; 
and  second,  a  consideration  of  the  combined  effects  of  Mach  number  and 
gaseous  imperfections,  with  principal  emphasis  in  the  latter  regard 
placed  on  the  caloric  imperfections  previously  discussed. 


Ideal  Gas  Flows 


Ihe  effects  of  Mach  number  of  primary  interest  here  are,  of  course, 
those  which  result  from  interaction  between  the  leading -edge  (or  other) 
shock  wave  and  small  disturbances  originating  on  the  surface  of  an  air¬ 
foil.  Some  insight  into  the  nature  and  extent  of  these  effects  can  be 
obtained  in  the  region  just  downstream  of  the  shock  wave  without  regard 


^For  a  given  value  of  T0,  Tpf,  to  the  accuracy  of  this  analysis,  is  the 
ideal  gas  function  of  Mq8n.  Thus,  knowing  %,  7n  can  be  determined. 

The  average  value  of  y  used  is  ya  =  7a(MQ5N)  =  ~  • 

6This  fact  was  employed  by  Linnell  (reference  15)  to  obtain  an  expression 
for  pressure  coefficient  equivalent  to  equation  (37)  for  the  case  of 
constant  7,  and  to  obtain  explicit  solutions  for  the  lift,  drag,  and 
pitching -moment  coefficients  of  several  airfoils  at  hypersonic  speeds. 
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for  the  shape  of  the  airfoil  producing  the  shock.  To  this  end  it  is 

bbfdCi  5p/SCi  .  »  %  N 

convenient  to  consider  the  ratio  - - - - -  (see  equation  (25)) 

b&/bc2  bp 

which  may  he  termed  "the  disturbance  strength  ratio"  since  in  the  region 
under  consideration  it  is  a  measure  of  the  ratio  of  strengths  of  distur¬ 
bances  reflected  from  the  shock  wave  to  disturbances  incident  on  the 
wave.  This  ratio  may  be  evaluated  with  the  expression 


(  gyNpN  *P  I  A 

db/dCi  _  Vsin  2  Pjj  d6  |  jj  /  sin  (%  +  “  °) 

ds/dc2  f  27NPN  ,  IP  I  \  sin  (%  -  5N  +  a) 

\sin  2  p,j  dS  |  N/ 


(38) 


which  is  easily  obtained  with  the  compatibility  equations  and  the 
oblique  shock-wave  equations  as  shown  in  appendix  D.  This  calculation 
has  been  carried  out  for  Mach  numbers  from  3*5  ho  »  (yN  =  1.4)  and 
flow  deflection  angles  approaching  those  corresponding  to  shock  detach¬ 
ment  (i.e.,  MN  %  l)  and  the  results  are  presented  in  figure  1.  It  is 
evident  that  except  near  Mjy  *  1,  the  ratio  is  small  (in  absolute  value) 
compared  to  1  throughout  the  entire  range  considered  -  this  observation 
also  applies,  of  course,  at  lower  supersonic  Mach  numbers*  Thus  it  is 
indicated  that  almost  all  of  an  incident  disturbance  is  generally- 
absorbed  in  the  shock  wave,  provided  the  air  behaves  like  an  ideal  dia¬ 
tomic  gas.7  This  result  is  substantially  the  same,  of  course,  as  that 
which  is  assumed  in  deriving  the  shock-expansion  method  of  calculating 
flows  about  airfoils,  and  therefore  yields  some  credence  in  the  method 
for  high  Mach  number  as  well  as  low  Mach  number  applications. 

As  an  over-all  check  on  the  shock-expansion  method,  surface  pres¬ 
sure  distributions  calculated  thereby  are  compared  in  figure  2  with 
those  obtained  with  the  method  of  characteristics  for  a  10 -percent -thick 
biconvex  airfoil  (a  =  0°)  operating  at  free-stream  Mach  numbers  of  3*5, 
5,  7*5,  10,  15,  and  <».  Predictions  of  the  slender-airfoil  approximation 
to  the  former  method  for  high  supersonic  speeds  are  also  shown.  There 
is  no  apparent  difference  between  the  ‘pressure  distributions  given  by 
the  method  of  characteristics  and  the  shock -expans ion  method  up  to  a 


7This  result  is  contrary  to  that  obtained  by  Ligh thill  (reference  l6) 
who  reports  that  for  hypersonic  flows,  a  disturbance  is  reflected 
from  a  shock  wave  with  opposite  sign  but  essentially  undiminished 
strength.  LighthillTs  conclusion  appears  to  be  based  on  an  incorrect 
evaluation  of  his  results  for  the  case  of  very  high  Mach  numbers . 
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Mach  number  of  10,  At  Mach  numbers  of  10,  and  above,  however,  the 
latter  method  predicts  pressures  which  are  slightly  low  downstream  of 
the  nose,  becoming  progressively  lower  with  increasing  Mach  number. 

This  result  would  be  deduced  from  figure  1  where  it  is  observed  that, 
at  the  Mach  numbers  under  consideration,  expansion  waves  incident  on 
the  nose  shock  are  reflected  back  toward  the  surface  as  compression 
waves  of  relatively  small  but  increasing  strength  with  increasing  Mach 
number.  The  effect  of  these  waves  does  not  become  pronounced  even  at 
infinite  Mach  number  (see  fig.  2(f))  and  the  shock -expans  ion  method  is 
thus  substantiated  as  being  a  reliable  simplified  method  for  predicting 
the  flow  about  airfoils  at  high  supersonic  speeds,  again,  so  long  as  the 
air  behaves  as  an  ideal  diatomic  gas.  The  further  simplified  slender- 
airfoil  method  also  appears  to  be  a  good  approximation  over  the  entire 
range  of  Mach  numbers,8  although,  as  would  be  expected  from  the  assump¬ 
tions  made  in  its  development,  it  is  in  somewhat  greater  error  than  the 
shock -expansion  method  at  the  lower  Mach  numbers. 

The  relative  accuracy  at  high  Mach  numbers  of  the  slender -airfoil 
method,  linear  and  second-order  potential  theories  may  be  seen  in 
figure  3*  As  might  be  expected,  the  slender-airfoil  method  is  more 
accurate  than  linear  theory  at  both  Mq  =  5  and  15*  and  more  accurate 
than  second -order  theory  at  Mq=15#  Itis  perhaps  surprising  to  note, 
however,  that  at  the  lower  Mach  number  of  5  the  slender-airfoil  method 
is  also  somewhat  superior  to  the  second-order  theory. 


The  pressure  distributions  of  figures  2  and  3  have  been  employed  to 
calculate  the  zero-lift  drag  of  the  biconvex  airfoil,  and  the  results  of 
these  calculations,  along  with  additional  predictions  of  linear  and 
second -order  theory,  are  shown  in  figure  k .  Predictions  of  the  shock- 
expansion  method  are,  of  course,  in  best  agreement  with  those  of  the 
method  of  characteristics;  while  the  slender-airfoil  method,  although 
slightly  less  accurate  than  the  shock -expans  ion  method,  is  apparently 
superior  to  both  linear  and  second -order  theory  at  Mach  numbers  above  3. 


The  preceding  findings  verify  that  so  long  as  the  disturbance 
strength  ratio  is  small  compared  to  1,  the  flow  along  streamlines  is 
essentially  of  the  Prandtl-Meyer  type.  If  we  choose,  on  the  basis  of 


these  findings,  a  maximum  absolute  value  for  of  0.06  (note 

^  58/5c2 


the  maximum  value  of 


56  /5c; 


for  the  cases  presented  in  fig.  2  was 


8The  hybrid  expression  for  pressure  coefficient  obtained  by  Ivey  and 
Cline  (reference  5)  gives  reasonably  good  results  also,  although  not 
as  accurate  as  the  slender-airfoil  method  at  the  higher  Mach  numbers 
under  consideration. 
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approximately  0.0 6  at  Mq  =  »),  the  region  in  which  the  shock-expansion 
method  is  applicable  can  readily  he  obtained  from  figure  1.  The  upper 
boundary  line  of  this  region  is  shown  in  figure  5  and  it  is  evident 
that  it  lies  only  slightly  below  (about  1°  in  general)  the  line  corre¬ 
sponding  to  shock  detachment  given  approximately  by  the  =  1.0  line. 
Almost  the  entire  region  of  completely  supersonic  (ideal  gas)  flow  is 
then  covered  by  the  method.9  (See  shaded  area  of  fig.  5.) 

The  question  naturally  arises  concerning  the  corresponding  range 
of  applicability  of  the  slender-airfoil  method.  This  question  may  be 
answered  in  part  by  comparing  separately  the  predictions  of  the  method 
for  oblique  shock  flows  and  expansion  flows  with  those  of  the  exact 
oblique  shock  equations  and  Prandtl -Meyer  equations.  Such  a  comparison 
is  shown  in  figure  6  in  terms  of  the  percentage  error  in  the  pressure 
coefficients  predicted  by  the  slender-airfoil  method.  As  would  be 
expected,  this  method  does  not  exhibit  good  accuracy  over  the  wide 
range  of  applicability  of  the  shock -expans  ion  method;  however,  it  is 
indicated  that  it  should  predict  pressure  coefficients  with  less  than 
10-percent  error  down  to  Mach  numbers  as  low  as  3  for, airfoils  producing 
flow  deflections  up  to  as  high  as  25  . 

As  a  further  check  on  the  utility  of  the  slender-airfoil  method, 
the  pressure  coefficients  on  the  10 -percent -thick  biconvex  airfoil  have 
been  calculated  with  this  method  and  the  shock-expansion  method  at  a 
Mach  number  of  10  and  angles  of  attack  up  to  about  30°#10  The  results 
of  this  calculation  are  shown  in  figure  7  (see  fig.  2(d)  for  a  =  0°) 
where  it  is  seen  that  the  agreement  is  reasonably  good  even  at  the 
highest  angle  of  attack.  This  fact  is  reflected  in  figure  8  showing 
the  force  and  moment  coefficients  for  the  airfoil  as  a  function  of  angle 
of  attack.  Little  difference  is  observed  in  the  force  coefficients  as 
calculated  by  the  two  methods,  while  the  moment  coefficients  display 
more  pronounced  but  nevertheless  small  differences  at  the  higher  angles 
of  attack. 


9If  it  is  required  as  by  Rand  (reference  17)  that  the  entire  flow  field 
be  of  the  true  Prandtl -Meyer  type  (i.e.,  that  all  flow  properties  be 
constant  along  first-family  Mach  lines  and  not  just  5  and  p) ,  then 
the  range  of  applicability  of  the  shock -expans  ion  method  would  be 
appreciably  smaller.  However,  it  has  been  shown  that  this  restriction 
is  not  necessary. 

10These  conditions  are  within  the  range  of  applicability  of  the  shock- 
expansion  method  as  defined  in  figure  4;  hence  the  use  of  the  method 
as  a  base  of  comparison  seems  justified.  Since  the  shock -expans ion 
method  is  far  less  tedious  to  apply  than  the  method  of  characteristics, 
it  will  be  employed  as  such  a  base  in  subsequent  calculations  whenever 
the  conditions  being  investigated  have  been  determined  to  be  within 
its  range  of  applicability. 
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From  these  and  previous  considerations,  the  ranges  of  applicability 
of  the  shock -expans ion  and  slender-airfoil  methods  for  ideal  supersonic 
gas  flows  are  reasonably  well  established.  It  remains  now  to  determine 
the  manner  and  extent  to  which  gaseous  imperfections  in  the  flow  at 

higher  supersonic  speeds  may  alter  these  ranges,  and  the  reasons 
therefor. 


Imperfect  Gas  Flows 


As  a  first  step  toward  investigating  the  effects  of  gaseous  imper¬ 
fections  on  the  high  Mach  number  flows  under  consideration,  it  is  con¬ 
venient  to  extend  our  consideration  of  the  disturbance  strength  ratio 
58/dCi 

—  .  It  is  recalled  that  when  air  exhibits  a  constant  value  of  y 

obfaC2 

equal  to  1.4  (the  value  for  an  ideal  diatomic  gas),  the  disturbance 
s  rength  ratio  is  small  at  arbitrarily  large  Mach  numbers,  provided  the 
flow  deflection  angles  are  not  too  close  to  those  for  shock  detachment. 
One  of  the  most  important  effects  of  gaseous  imperfections  is,  however 
to  decrease  y  of  the  disturbed  air  below  this  value  due  to  the  excita¬ 
tion  of  additional  degrees  of  freedom  (e.g.,  vibrational)  in  the  mole¬ 
cules  at  the  high  temperatures  encountered  at  high  Mach  numbers.  Indeed 
a  arbitrarily  high  Mach  numbers  it  might  be  expected  that  y  of  the 
isturbed  air  would  approach  1,  since  the  number  of  degrees  of  freedom 
may  effectively  become  very  large  (see,  e.g.,  references  3  and  6).  In 
this  case,  however,  the  extent  of  the  disturbance  flow  field  is  decreased 

+°  aLth€  surface  of  the  tody  which  is  negligibly  thin  compared 

to  that  for  the  case  of  ideal  gas  flow.  Thus  it  is  apparent  that  signi- 
can  changes  in  the  flow  about  airfoils  at  high  Mach  numbers  may 
result  from  decreases  in  y  of  the  disturbed  air;  hence  the  effects  of 

such  decreases  on  the  disturbance  strength  ratio  would  appear  to  warrant 
attention. 


A  detailed  analysis  of  these  effects  is  impractical  at  the  present 
time  due  to  the  limited  range  over  which  the  variation  of  y  with  tem¬ 
perature  is  known.  However,  some  knowledge  of  these  effects  can  be 
gained  by  repeating  the  ideal  gas  calculations  for  constant  values  of 
7jj  between  1.4  and  1.0. 11  Such  calculations  have  been  carried  out  at 


Since  the  enthalpy  is  negligibly  small  compared  to  the  mass  kinetic 
energy  of  the  undisturbed  fluid  at  the  high  Mach  numbers  of  interest 
and,  hence,  y  of  this  fluid  does  not  influence  the  flow,  this  approach 
corresponds  to  employing  an  average  value  of  y  for  the  disturbed 
fluid.  Since  only  flows  of  dense  air  are  considered  here,  heat- 
capacity-lag  phenomena  are  neglected  (see  references  5  and  6). 
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ds/dCi 

infinite  Mach  number  since  in  this  case  -  has  its  maxim-urn  value 

,  5s/dC2 

for  a  given  7^,  and  the  results  are  presented  in  figure  9*  It  is  seen 
that  except  near  shock  detachment,  the  disturbance  strength  ratio 
increases  with  decreasing  7^,  approaching  1  as  7^  approaches  1.  This 

increase  is  slow  at  first;  for  example,  the  value  of  is  still 

Sb  fdC2 

less  than  0.1  at  7^  =  1.3*  It  might  therefore  be  expected  that  the 
shock -expans ion  method  would  continue  to  apply  with  reasonable  accuracy 
so  long  as  7  of  the  disturbed  flow  is  not  appreciably  less  than  this 
value.  This  point  has  been  checked  with  the  methods  developed  pre¬ 
viously  for  analyzing  the  flow  of  a  calorically  imperfect  diatomic  gas 
at  local  air  temperatures  up  to  about  5000°  R  (note  7  has  a  value 
only  slightly  less  than  1.3  at  this  temperature).  In  particular,  the 
pressure  distribution  on  the  lower  surface  of  the  biconvex  airfoil  at 
M0  =  10,  a=19.9°,  and  T0  =  500°  R  (1^4000°  R  at  leading  edge)  has 
been  calculated  with  both  the  method  of  characteristics  and  the  shock- 
expansion  method.12  The  results  of  these  calculations  are  presented  in 
figure  10  and  it  would  appear  that  the  conclusions  drawn  from  figure  9 
pertaining  to  cases  where  7^  is  of  the  order  of  1.3  or  greater  are 
substantiated.  Pressures  in  the  expansion  flow  about  the  upper  surface 
are  not  influenced  (due  to  the  low  temperatures)  by  caloric  imperfections 
and  hence  are  the  same  as  shown  in  figure  7(b). 

"  Shown  also  in  figure  10  is  the  pressure  distribution  obtained  by 

the  shock -expans ion  method  for  an  ideal  gas  (7^  =  1.4).  It  is  apparent, 

W  on  comparing  this  pressure  distribution  with  the  other  distributions, 

that  although  the  effect  of  caloric  imperfections  on  the  disturbance 
strength  ratio  is  small,  the  pressures  are  appreciably  reduced  by  the 
increase  in  specific  heats.  The  extent  of  this  reduction  is  more  com¬ 
pletely  illustrated  in  figure  11  where  the  lower  surface  pressure  distri¬ 
butions  on  the  biconvex  airfoil  are  presented  for  Mo=10  and  To=500°  R, 
at  a  =  0°,  10°,  19.9°,^  and  30°.  As  one  might  expect,  the  reduction  in 
pressures  increases  with  angle  of  attack  (due  to  the  corresponding 
increase  in  static  temperature  of  the  disturbed  air) •  The  pressure 
coefficients  calculated  with  consideration  for  the  imperfections  in  the 
gas  are  less  on  the  lower  surface  (up  to  6  percent  at  the  leading  edge 
and  15  percent  at  the  trailing  edge)  than  those  calculated  assuming  the 
gas  behaves  ideally.  The  upper-surface  pressures  are  again  unaffected 

12For  added  ease  of  calculation  the  expansion  method  of  appendix  C  was 
employed.  This  method  is  also  employed  in  all  subsequent  calculations 
of  this  type  since  it  has  been  found  to  yield  results  differing  by 
s  less  than  1  percent  from  those  obtained  by  the  more  tedious  graphical 

integration  method. 
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by  the  caloric  imperfections  of  air  in  all  the  cases  presented  (except 
at  a  -  0  )  since  this  surface  experiences  lower  pressures  and  hence 
lower  temperatures.  They  are  therefore  the  same  as  shown  in  figure  7. 
Shown  also  in  figure  11  are  the  pressure  distributions  calculated  with 
the  slender-airfoil  method  for  y  =  7a.  The  accuracy  of  this  simpli¬ 
fied  method  is  substantially  the  same  as  was  previously  observed  for 
the  corresponding  method  in  the  case  of  ideal  gas  flows,  although  the 
local  error  may  be  greater  than  the  reduction  in  pressure  coefficients 
due  to  the  caloric  imperfections  of  air.  This  error  is  somewhat  com¬ 
pensating,  however,  in  its  effects  on  the  force  and  moment  coefficients, 
as  will  be  seen.  * 

The  force  and  moment  coefficients,  corresponding  to  the  lower- 
surface  pressure  distributions  shown  in  figure  11  and  the  upper-surface 
distributions  of  figure  7  are  presented  in  figure  12.  The  reduction  in 
the  lower-surface  pressures  leads,  of  course,  to  a  general  reduction  in 
all  three  coefficients  (up  to  about  10  percent  for  a  =  30°) .  The 

slender-airfoil  method  again  predicts  these  coefficients  with  surprising 
accuracy.  *  ° 

In  order  to  further  assess  the  accuracy  of  the  slender-airfoil 
method  some  additional  calculations  were  carried  out  for  the  biconvex 
airfoil  at  a  =  0°  and  Mq  =  20  and  30.  The  pressure  distributions  for 
these  cases  were  calculated  by  the  shock-expansion  method,  slender- 
airfoil  method  (7  =  7a),  and  slender-airfoil  method  (7  =  7^).  These 
results  are  presented  in  figure  13  and  it  is  observed  that  the  use  of 
7a  rather  than  y^  improves  the  accuracy  of  the  slender-airfoil  method. 
The  extent  of  this  improvement  in  the  case  of  drag  coefficient  is  shown 
in  figure  14  “it  would  appear  that  predictions  of  the  slender-airfoil 
method  (7  =  7a)  and  shock-expansion  method  are  in  as  good  agreement  as 
for  ideal  gas  flows  (see  fig.  4).  On  the  basis  of  these  and  previous 
results,  it  may  be  concluded  then  that  not  only  does  the  shock-expansion 
method  retain  its  range  of  applicability  when  air  exhibits  caloric 
imperfections  provided  7  of  the  disturbed  air  is  not  appreciably  less 
than  1.3,  but  also  the  slender-airfoil  method  (7  =  ya )  retains  its 
range.  “* 

It  would  be  surprising  indeed,  however,  if  this  conclusion  con¬ 
tinued  to  apply  as  7  of  the  disturbed  fluid  approached  1  since,  as 

discussed  previously,  —  ■■  1  is  not  small  compared  to  1  in  this  case, 

36  /dc2  1 

but  would  appear,  in  fact,  to  approach  1.  This  matter  may  be  investi¬ 
gated  in  the  same  manner  as  the  effect  of  7N  on  the  disturbance 

strength  ratio  was  investigated,  namely,  by  using  the  ideal  gas  relation¬ 
ships  in  combination  with  appropriate  values  of  7. 
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The  limiting  case  of  infinite  free -stream  Mach  number  and  7  =  1.0 
(for  the  disturbed  fluid,  see  footnote  11)  has  already  been  investigated 
by  Busemann  (reference  18)  and  more  recently  by  Ivey,  Klunker,  and  Bowen 
(reference  19)*  In  this  case,  as  pointed  out  previously,  the  shock  wave 
emanating  from  the  leading  edge  remains  attached  to  the  surface  down¬ 
stream  of  the  leading  edge  (this  is  easily  verified  with  the  oblique 
shock-wave  equations)  and  the  disturbance  flow  field  is  confined  to  an 
infinitesimally  thin  layer  adjacent  to  the  surface.  In  addition,  the 
velocity  along  a  streamline  downstream  of  the  shock  is  constant,  as  may 
easily  be  shown  with  the  compatibility  equations.  Surface  pressures 
therefore  become  a  simple t function  of  airfoil  geometry, 

d6q 

Cp  =  2  sin2  6g  +  2  cos  5g 


varying,  to  a  first  approximation,  directly  with  the  square  of  the  com¬ 
ponent  of  free-stream  velocity  normal  to  the  surface  (i.e.,  the  flow  is 
approximately  of  the  Newtonian  corpuscular  type).  With  this  theory 
then,  and  the  method  of  characteristics,  we  can  get  an  idea  of  both  the 
extent  to  which  extreme  changes  of  7  from  1.4  toward  1  will  alter 
surface  pressures,  and  the  accuracy  with  which  the  shock -expans ion 
theory  predicts  the  alterations.  To  this  end,  figure  15  is  presented 
showing  the  pressure  distributions  about  the  biconvex  airfoil  at 
Mq  =  00  as  calculated  by  the  several  methods  for  different  values  of  7* 
It  is  observed  that,  whereas  the  shock -expans ion  method  agrees  very 
closely  with  the  method  of  characteristics  for  7  =  1.4,  there  is  a 
large  difference  at  7  =  1.05-  This,  of  course,  is  precisely  what  one 
would  expect  from  the  previous  discussion  of  the  disturbance  strength 
ratio.  On  the  other  hand,  if  the  two  characteristic  solutions  and  the 
Busemann  method  are  considered  in  order  of  decreasing  7 >  it  is  indi¬ 
cated  that  the  characteristics  solutions  approach  the  Busemann  theory 
as  7  approaches  1.  For  7  =  1.0  and  M 0  =  00  the  shock-expansion 
method,  in  turn,  predicts  a  discontinuous  pressure  distribution  with  a 
pressure  coefficient  equal  to  that  of  the  Busemann  theory  at  the  leading 
edge  but  a  pressure  coefficient  of  zero  at  all  points  downstream  of  the 
leading  edge.  Hence  it  may  be  concluded  that  when  the  free-stream  Mach 
number  approaches  infinity  and  7  approaches  1,  the  Busemann  method 
rather  than  the  shock -expans  ion  method  for  calculating  the  flow  about 
airfoils  should  be  employed. 

CONCLUDING  REMARKS 

The  flow  about  curved  airfoils  was  investigated  analytically  at 
high  supersonic  speeds  first  assuming  air  behaves  as  an  ideal  gas,  and 
then  assuming  air  behaves  as  a  thermally,  although  not  necessarily 
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calorically,  perfect  gas.  (Caloric  imperfections  had  previously  been 
observed  to  be  of  predominant  importance  in  free  flight,  at  least  for 
local  air  temperatures  up  to  the  order  of  5000°  R.)  It  was  found  that  « 

so  long  as  air  exhibits  no  imperfections  (i.e.,  behaves  as  an  ideal 
diatomic  gas)  small  disturbances  originating  on  the  curved  surface  of 
an  airfoil  are  almost  completely  absorbed  in  the  shock  wave  emanating 
from  the  leading  edge  (or  any  other  location  on  the  surface),  provided 
he  flow  deflection  angles  are  not  too  close  to  those  corresponding  to 
shock  detachment.  This  result  is  essentially  independent  of  Mach  number 
and  is  consistent  with  the  early  calculations  of  Crocco  concerning  sur-  * 
face  pressure  gradient  at  the  leading  edge  of  an  airfoil.  It  was  con¬ 
cluded  that  in  ideal  gas  flows  the  shock -expans ion  method  for  determining 
flow  conditions  at  the  surface  of  an  airfoil  would  apply  with  good 
accuracy  at  arbitrarily  high  Mach  numbers.  This  conclusion  was  verified 
by  the  excellent  agreement  found  between  pressure  distributions  on  a 
10 -percent -thick  biconvex  airfoil  at  0°  angle  of  attack  calculated  with 
the  shock -expans ion  method  and  the  method  of  characteristics.  It  was 
further  shown  that  the  former  method  can  he  easily  employed  to  construct 
e  entire  flow  field  about  an  airfoil  in  a  manner  that  would  account 
for  shock-wave  curvature  and  entropy  gradients  resulting  therefrom. 


A.higl?  Mach  n™ber  ^  approximation  to  the  shock -expans  ion  method  for 
in  airfoils  was  also  investigated,  and  was  found  to  apply  with  good 
accuracy  at  Mach  numbers  above  3  and  flow  deflection  angles  up  to  25°. 
The  essential  feature  of  this  slender-airfoil  method  is,  of  course,  its 

simplicity,  and  for  that  reason  it  may  prove  useful  for  some  engineering 
purposes. 


Effects  of  caloric  imperfections  were  first  investigated  aualita- 

^TtLrr  nrirVhe  reaU0tlOn  “  the  ratl°  specifi?  hefts  from 

1.4  toward  1  0.  It  was  found  that  as  the  ratio  decreased,  the  extent 

SmS  «  ?h1S  ”)“CeS  refleCtel  fr0m  a  sh“k  —  ^creased  iflL 

limit  as  the  ratio  approached  1,  the  reflection  was  conplete,  and  the 
raff  IT  TTS  tangent  to  the  surface  of  the  airfoil.  So  long  as  t 

less°than  lo below  1.3,  however,  it  was  foLl  th 
10  percent  of  a  disturbance  was  reflected;  hence  the  simple 
shock-expansion  method  might  be  expected  to  continu^  to  apply.  This 
matter  was  checked  quantitatively  for  the  biconvex  airfoil  with  the  ai 

imDeT-?611^112^  m®thod  of  characteristics  including  effects  of  calori 
imperfections  (up  to  local  air  temperatures  of  the  order  of  5000°  R 

°  a  ratl°  0f  k“ts  of  about  1.3).  Itwas  t’o  un, 

irlaawL  T ck-expansion  method  was  substantially  as  accurate  as  for 
JT  g  provided  lb  also  generalized  to  include  effects  o' 

heat  raSor^  t°nS'^  principal  effect  of  the  reduction  in  specif i, 
gas  valued  r,!  T*  T  pressure  coefficients  below  their  ideal 
coefficier.TJ!’  much  as  15  percent.  The  reduction  in  force  and  momen- 
it  1  ™  was  somewhat  smaller,  being  about  10  percent.  Similarly 

it  was  found  that  the  slender-airfoil  method,  modified  to  employ  an 
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average  value  of  the  ratio  of  specific  heats  for  a  particular  flow 
field,  exhibited  essentially  the  same  accuracy  as  the  analogous  method 
for  ideal  gas  flows.  Thus  it  is  indicated  that  the  generalized  shock- 
expansion  method  and  its  slender-airfoil  counterpart  can  he  applied 
with  good  accuracy  up  to  very  high  free -flight  Mach  numbers.  If  the 
flow  deflection  angles  are  less  than  about  12°,  these  Mach  numbers  are 
of  the  order  of  the  so-called  escape  Mach  number  (i.e.,  Mach  numbers 
as  high  as  30  to  35) . 

At  even  higher  Mach  numbers  where  the  ratio  of  specific  heats  is 
expected  to  decrease  appreciably  below  1.3,  and  in  fact  perhaps  to 
approach  1,  it  was  not  possible  to  obtain  an  accurate  check  on  the  pre¬ 
viously  discussed  qualitative  considerations.  It  was  undertaken,  how¬ 
ever,  to  compare  the  pressure  distributions  on  the  biconvex  airfoil 
predicted  by  the  Busemann  method  (for  the  limit  of  the  ratio  of  specific 
heats  approaching  1  and  Mach  number  approaching  infinity)  with  those 
predicted  by  the  shock-expansion  method  and  the  method  of  characteris¬ 
tics  at  infinite  Mach  number  (employing  a  constant  ratio  of  specific 
heats  of  1.05)  •  As  was  expected,  the  shock -expans  ion  method  was  in 
very  poor  agreement  with  the  method  of  characteristics,  whereas  the~ 
Busemann  method  was  in  relatively  good  agreement.  It  is  therefore 
indicated  that  for  extremely  high  Mach  numbers  (something  in  excess  of 
the  sea-level  escape  Mach  number)  the  Busemann  theory  may  apply. 
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APPENDIX  A 


METHOD  OF  CHARACTERISTICS  FOR  TWO-DIMENSIONAL  FLOW 
OF  A  CALORICALLY  IMPERFECT  GAS 


In  the  application  of  the  method  of  characteristics  for  a  calori- 
cally  imperfect  diatomic  gas  to  the  particular  problem  of  analyzing  the 
flow  about  curved  two-dimensional  airfoils,  many  of  the  calculations 
are  identical  to  those  encountered  in  the  solution  of  any  problem  where 
characteristics  theory  is  employed.  Since  the  details  of  these  calcu¬ 
lations  are  well  known  and  well  reported  (see,  e.g.,  reference  ll), 
they  will  not  be  repeated  here. 

A  lattice -point  system  with  an  initial -value,  numerical  computing 
procedure  will  be  used.  The  form  of  the  compatibility  equations  to  be 
employed  was  developed  previously;13  however,  it  is  convenient  for  pur¬ 
poses  of  calculation  to  substitute  the  pressure  ratio,  p/q^  into  these 
equations  and  to  rewrite  them  as  difference  equations.  Equations  (12) 
and  (13)  are  thus  reduced  to  the  following  forms 

(p/qG)c  -  (p/q0)A  =  -Xa(5C  -  5A)  (Al) 

and 

(p/qo)c  -  (p/qo)B  =  (5c  -  &b)  (A2) 

where 

x  = 

sin  23  (A3) 

It  is  also  convenient  to  employ  several  reference  curves.  These  curves 
can  he  divided  into  two  groups.  The  general  reference  curves  consist  of 
y  and  t(T)  as  a  function  of  temperature,  T.  Equations  (l4)  and  (17)  are 

Q  '  i'i  —  . 

This  form  of  the  compatibility  equations  (in  p  and  8  coordinates) 
was  also  used  in  obtaining  some  of  the  characteristics  solutions  for 
ideal  gas  flows.  The  majority  of  these  solutions  were  carried  out, 
however,  with  the  compatibility  equations  in  (3,  8,  and  entropy 
coordinates,  since  it  was  found  that  greater  accuracy  was  usually 
obtained  for  a  given  net  size.  In  general,  the  net  size  employed 
yielded  pressures  at  from  30  to  35  surface  points  on  an  airfoil  with 
a  maximum  error  in  the  corresponding  pressure  coefficients  equal  to 
less  than  1  percent  of  the  pressure  coefficient  at  the  leading  edge. 
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used  to  determine  these  curves.  A  second  set  of  shock -'wave  reference 
curves  consisting  of  p/qQ,  and  &  as  a  function  of  temperature,  T, 
are  determined  by  use  of  equations  (l8)  through  (20 )  -  the  values  of  TQ 
and  Mo  are  presumed  known. 

In  the  computations  three  types  of  points  are  encountered.  These 
are  (l)  a  point  in  the  flow  field  between  the  shock  wave  and  the  air¬ 
foil  surface,  (2)  a  point  on  the  airfoil  surface,  and  (3)  a  point  just 
downstream  of  the  shock  wave.  Each  one  of  these  types  of  points 
requires  a  slightly  different  computing  procedure  and  they  will  be  con¬ 
sidered  in  order. 


Point  in  the  Flow  Field  Between  the  Shock  Wave 
and  the  Airfoil  Surface 

i 

Figure  16(a)  shows  a  schematic  diagram  of  the  system  of  points  to 
be  considered  in  these  calculations.  Point  C  is  the  unknown  point  at^_ 
the  intersection  of  the  first-family  characteristic  line  passing  through 
point  A  and  the  second -family  characteristic  line  passing  through 
point  B.  Six  quantities  are  known  at  both  points  A  and-  B,  and  the 
problem  is  to  calculate  these  same  quantities  at  point  C.  These  quanti¬ 
ties  are  x,  y,  5,  p/q^  T,  and  %.  The  first  five  quantities  are  of 
obvious  significance.  The  sixth,  Tyj,  is  defined  as  the  static  tempera¬ 
ture,  just  downstream  of  the  shock  wave,  on  the  streamline  passing 
through  the  point  C. 


The  physical  coordinates  of  the  point  C  (xq,  yp)  may  be  determined 
by  standard  procedures  such  as  those  given  in  reference  11.  In  order  to 
determine  the  quantity  5q,  it  is  necessary  to  solve  equations  (Al)  and 
(A2)  simultaneously;  thus 


_  XA  5  A  +  %  5B  +  fo/Op )A  ~ 


XA  +  \B 


(A4) 


Equation  (Al)  or  (A2)  is  then  used  to  obtain  (p/q,-,)^  • 

There  remains  only  the  problem  of  determining  Tq  and  T™  at 

point  C.  The  temperature  T^  is  obviously  constant  along  the  stream- 

C 

line  through  C.  This  quantity  may  therefore  be  calculated  in  the  same 
manner  as  the  entropy  is  calculated  in  similar  flow  fields  for  ideal  gas 
processes  (see,  e.g.,  reference  11).  Furthermore,  since  the  flow  along 
streamlines  downstream  of  the  shock  wave  is  isentropic,  equation  (l6) 
may  be  applied  in  the  form 
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(p/9o)c  ^(Tnc) 

(p/<1o)n  ’  *(TC) 

0 


The  pressure,  (p/q^  ,  is  defined  in  a  manner  analogous  to  TNc  ,  and 

may  thus  he  determined  using  the  shock -wave  reference  curves  and  the 
known  value  of  •  Similarly,  ^  )  may  he  determined  from  the 


general  reference  curves.  The  only  unknown  in  equation  (a5)  then  is 
y  (Tc)  which  may  now  he  calculated.  Once  (Tq)  is  determined,  Tq  may 
he  determined  hy  again  using  the  general  reference  curves.  All  six 
quantities,  xQ,  yc,  5C,  (p/go)0  ,  Tc,  and  TN  have  now  heen  determined. 


Point  on  the  Airfoil  Surface 


Figure  16(b)  shows  a  schematic  diagram  of  the  points  to  he  con¬ 
sidered  in  these  calculations.  The  physical  coordinates  of  point  C, 
(XC>  7c)*  are  firs't  calculated  hy  solving  simultaneously  the  equation 
of  the  second -family  Mach  line  passing  through  point  B  and  the  equation 
of  the  airfoil  surface.  When  xc  and  yc  have  heen  determined,  8r  is 


readily  obtained  from  the  equation  of  the  airfoil  surface, 
is  then  applied  to  determine  (p/qQ) 


Equation  (A2) 


Since  the  airfoil  surface  is  a  streamline, T^  is  constant  along  the 

surface  and  may  he  evaluated  at  the  leading  edge.  The  temperature,  Tq, 

may  then  he  determined  using  equation  (A5)  and  the  previously  described 
procedure.  All  six  quantities,  xc,  yc,  (p/qQ)  ,  8C,  Tc,  and  ,  are 
thus  determined.  L  L 


In  the  special  case  of  the  first  point  on  the  airfoil  surface 
downstream  of  the  leading  edge,  the  pressure  ratio  is  calculated  using 
the  procedure  of  reference  9«  This  procedure  is  easily  shown  to  he 
applicable  to  calorically  imperfect  gas  flows  providing  the  oblique 
shock -wave  equations  of  the  present  paper  are  employed. 


Point  on  the  Shock  Wave 


Figure  16(c)  shows  a  schematic  diagram  of  the  points  to  be  con¬ 
sidered  in  these  calculations.  The  physical  coordinates  of  point  C 
(XC>  yc)  are  first  calculated  hy  solving  simultaneously  the  equation 
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of  the  first-family  Mach  line  passing  through  point  A  and  the  equation 
of  the  shock  wave  linearized  at  point  D,  the  last  known  point  on  the 
wave.  The  variation  of  p/qQ  with  5  along  the  shock  wave  may  he 
approximated  by  the  relation 


(P/K>>c  '  ~ 


In  this  equation 


C 

*(p/gp) 

as 


(Sf!  “  &I)) 


N 


(AS) 


is  the  rate  of  change  of  p/qQ  with  5 


N 


along  the  downstream  side  of  the  shock  wave  evaluated  at  point  D. 
Because  of  the  complicated  nature  of  the  shock -wave  equations,  it  is 

d(p/q.0) 

generally  easier  to  evaluate 


dS 


graphically  or  numerically 


N 


from  the  shock -wave  reference  curves.  Equations  (Al)  and  (A6)  are 
solved  simultaneously  for  Sq  ,  thus 


a(p/qQ) 


5C  = 


XA  5A  +  d5 


“  '  8D  +  (p/<Jo)a  "  (pALo)d 


i(p/Qo) 

dS 


(A7) 


N 


When  Sc  has  been  calculated,  Tq  ,  and  in  turn  (p/qo)^  ,  may  be  deter- 

mined  from  the  shock -wave  reference  curves.  Since  point  C  in  this  case 

is  just  downstream  of  the  shock  wave,  Tc  and  are  identical. 

0 

The  six  quantities ,  Xq,  y q,  (p/<l0)  *  9  an^  TN n  ^ave  now  ^een 

C  C 

determined. 
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APPENDIX  B 

SHOCK-EXPANSION  METHOD  FOR  CALCULATING  THE  FLOW  * 

FIELD  ABOUT  AN  AIRFOIL 


An  initial  value  procedure  which  is  similar  to,  although  markedly 
simpler  than,  that  associated  with  the  method  of  characteristics  may 
he  employed  to  carry  out  this  calculation.14-  To  illustrate,  consider 
the  sketch:  9 


With  the  oblique  shock -wave  and  expansion  equations,  all  fluid  proper¬ 
ties  at  points  M,  A,  C,.  and  so  forth  on  the  airfoil  surface  may  he  cal¬ 
culated  in  the  usual  manner.  If  the  point  A  is  chosen  close  to  M,  the 
first-family  (Cx)  Mach  line  connecting  A  to  point  B  on  the  shock  wave 
may  he  considered  straight  and  inclined  at  an  angle  to  the  free -stream 
direction  equal  to  p^  4-  5^.  Similarly,  the  segment  MB  of  the  shock 
wave  may  he  considered  straight  and  inclined  at  the  angle  crM  to  the 
stream  direction.  Thus  the  point  B  in  the  flow  field/ may  easily  he 
determined.  The  direction  of  BD  (a  segment  of  the  streamline  passing 
through  B)  is  the  same  as  the  tangent  to  the  surface  at  A,  and  the 
attitude  of  the  segment  BE  of  the  shock  wave  is  fixed  hy  this  direction* 
The  locations  of  points  D  and  E  in  the  flow  field  are  thus  fixed  once 
point  C  is  chosen.15  The  construction  of  the  remainder  of  the  flow  field 
follows  in  a  similar  fashion.  Having  determined  the  shapes  of  the 
streamlines,  the  fluid  properties  along  these  lines  are,  of  course, 
determined  in  the  same  manner  as  those  along  the  surface. 


14It  is  clear,  of  course,  that  an  ”average  value”  procedure  could  also 
he  employed.  Such  a  procedure  would,  in  fact,  he  the  more  desirable 
l5  some  cases,  since  a  coarser  net  may  he  used. 

The  point  C  should,  of  course,  he  chosen  close  to  A. 
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It  is  important  to  note  that  this  method  is  applicable  to  the 
determination  of  the  flow  not  only  in  the  region  adjacent  to  the  air¬ 
foil  (whether  the  surface  he  concave  or  convex)  hut  also  in  the  region 
downstream  of  the  airfoilj  hence  it  may,  for  example,  prove  useful  in 
downwash  studies  and  the  like. 
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APPENDIX  C 

APPROXIMATE  SOLUTION  FOR  PRANDTL -MEYER  FLOW 
OF  A  CALORICALLY  IMPERFECT  GAS 


The  following  solution  is  obtained  with  an  analysis  similar  to 
that  used  in  Meyer’s  original  paper  (reference  20).  A  schematic  dia¬ 
gram  of  the  subject  flow  field  is  shown  in  figure  17.  It  is  evident 
that  the  change  in  flow-inclination  angle  for  Prandtl -Meyer  flow  can 
be  written  as  follows 


6N  "  8  =  (%  "  0)  +  (WR  "  w) 


(Cl) 


Since  the  flow  is  isentropic,  a  given  value  of  the  local  pressure  will 
determine  the  Mach  angle,  p.  The  problem  then  is  to  evaluate  the 

To  this  end  the  velocity  components  tangential  and  normal 
to  the  first-family  Mach  lines  may  be  expressed  in  the  usual  manner  in 
terms  of  a  potential  cp,  thus 


dtp 

Sr 


(C2) 


_  1  ST 

r  Sgj 


(C3) 


It  is  clear,  however,  that  these  components  are  functions  of  w  only; 
hence  it  is  convenient  to  define  a  new  velocity  potential  which  is  a 
function  of  u  alone.  Such  a  potential  is 


(c4) 


The  velocity  components  may  then  be  written  in  terms  of  this  new 
potential. 


w 


u  =  $ 

(C5) 

a  = 

(c  6) 

The  resultant  velocity  is  given  by  the  expression 

v2  =  o2  +  qJ  (C7) 
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Equation  (15)  for  conservation  of  energy  may  be  written  in  terms  of  the 
local  temperature  as  follows: 


V2  +  2  (  JRT  +  2RT 


Vee/T_i 


The  constant.  A,  is  evaluated  at  the  conditions  existing  upstream  of  the 
expansion  region;  namely > 


A2  =  Vfl2  + 


+  2RTN 


e/% 


eD/  iK_i 


Equations  (C7)  and  (C8)  are  then  combined  to  yield 


$2  +  =  -2RT  (22l_  +  — £/2L  ^  +  A2 

Vn-1  e0/T-]V 


(CIO) 


It  was  shown  previously,  however,  that 


a2  =  7RT 


(oil) 


Equations  (C6),  (CIO),  and  (Cll)  may  therefore  be  combined  to  obtain 
the  following  relationship: 


ee/T_x 


(C12) 


+  ,  a  (Zili+  TZi  .!  +  Zi±  (Z±  W*_ 
.  “  l  7i-l  7±-l  L  7  7i  \  7  y  e0/T- 1 


=  A2  (C13) 


From  the  imperfect  gas  relationship  for  7  we  have 


1  +  (V1}  \  T  /  (e0/T-l)2 


(e0/T-l)2 


(Cl4) 


1 
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Substituting  this  relation  into  equation  (C13)  there  is  then  obtained 


$2 


+ 


0) 


(C15) 


* 


where 


F 


i)- 


e/T 


1  +  (74-1) 


(e/T)- 


e/T 


(ee/T-l) 


(ee/T-l) 


1+ 


6/T 


(e^-Dj 


1  + 


(71-1)  (e/T)2e0/T 
7i  (e0/T-l)2 


(ci6) 


Now 


7T  =  a2 

7*0  7^0 


(C17) 


For  every  value  of  T/0  there  is  thus  a  particular  value  of  as/71R0. 
The  function  F ( 0/T )  is  therefore  uniquely  determined  for  any  value 

of  a2  since  7*R0  is  of  course  a  constant.  With  this  point  in  mind, 
let  ' 


F(0/T)  =  G(aa/71R0)  (Cl8) 

Figure  18  shows  G(a2/7*R0)  plotted  as  a  function  of  a2/7*R0.  This 
curve  is  approximated  with  the  following  simple  relation: 


G(a2/7iR0) 


for  0.18  < 

7±R& 

and 


<  1.0 


°*38  +  0.71 

7jR0 


0.l4 

a2/ 7*R0 


a2 

7iR0 


G(a2/71R0)  =  0 


(ciy) 


for  0  < 


<  0.18 


(C20) 
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Equation  (C23)  then  reduces  to 

4  /  \4  4  COS  v  O 

COS4  T  ( T(jj) 4  +  - -  COS1 2 3 4 5 6 7  T  (TJ2  + 

D 

This  equation  is  solved  for  tu,  thus 


4  sin2  t 
D2 


4  sin2  v 


(C28) 


tr 


Tw  =  '  - 


2  1 


D  cos  T 


(cos  T  -  COS  V 


(C29) 


or 


du 


■A 


COS  T  dT 


(cos  T  -  cos  v)1/2 


(c  30) 


This  egression  is  readily  integrated  to  obtain  the  following  equation 
relating  m  to  the  local  velocity: 


to  -  Wjj  =  D  -^2 


F(k,z)  —  E(k,z^) 


F(k,z) 


where 

E  elliptic  integral  of*  the  second  kind 
F  elliptic  integral  of  the  first  kind 
k  sin  ~  (modulus) 

z  sin-1  ('-l1?  )  (amplitude) 

4  sin  v/2  / 


F(k,  zN) 


1 

f 


(C3l) 


The  procedure  for  calculating  corresponding  values  of  the  pres 
sure  p,  and  the  deflection  angle,  8,  is  straightforward  with  the  a 
he  preceding  equations  and  may  be  summarized  as  follows: 


1.  Calculate  A2,  equation  (C9) 

2.  Calculate  D2,  equation  (C24) 

3.  Calculate  v,  equation  (C25) 

4.  Assume  a  value  of  T,  less  than  T„ 

5.  Calculate  p,  equations  (l6)  and  (17) 

6.  Calculate  V2,  equation  (C8) 

7.  Calculate  y ,  equation  (Cl4) 
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8.  Calculate  a2  (or  S^2) ,  equation  (Cll) 

9.  Calculate  M  and,  in  turn,  3  from  V  and  a 

10.  Calculate  u2  (or  d>2) ,  equation  (C7) 

11.  Calculate  t,  equation  (C26) 

12.  Calculate  w  ,  equation  (C3l) 

13.  Calculate  8,  equation  (Cl) 

it. is  procedure  is  followed  so  long  as  the  quantity  a2/ y ^R0  is 
greater  than  0.l8.  (This  is  equivalent  to  the  temperature  being  greater 
than  approximately  1000°  R.)  For  values  of  a2/y1R0  less  than  or  equal 
to  0.18  (or  temperatures  less  than  about  1000°  R)  G  is  set  equal  to  zero 
(see  equation  (C20)).  In  this  case  equation  (C15)  reduces  to  the  same 
form  as  for  an  ideal  gas,  and  therefore  the  well-known  ideal  gas  rela¬ 
tionships  can  be  used. 
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APPENDIX  D 

EVALUATION  OF  THE  DISTURBANCE  STRENGTH  RATIO 

Consider  the  element  of  an  oblique  shock  wave  shown  in  the  accom¬ 
panying  sketch 


strengths  of  reflected  and  incident  disturbances,  respectively,  since 
all  disturbances  incident  on  the  wave  between  points, D  and  C  must 
travel  along  first -family  characteristics  which  cross  C2  between 
D  and  A,  while  all  disturbances  reflected  from  the  shock  wave  in  this 
region  must  travel  along  second-family  characteristics  which  cross  Cj_ 
between  A  and  C.  This  ratio,  termed  the  disturbance  strength  ratio, 
may  be  evaluated  locally  just  downstream  of  the  shock  wave  in  the  fol¬ 
lowing  manner.  The  points  D  and  C  are  chosen  sufficiently  close 
together  so  that  the  difference  in  pressure  between  these  points  may 
be  written 


or 


Pc  -  Pp 


?  6p 

r~  ac2  +  c~~  ac 
a  6cx 


Pc  "  Pd 


dp 

—  (5c-Sd) 

dS  N 


(PI) 

(D2) 
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where  t?  is  the  rate  of  change  of  pressure  with  flow  deflection 
d6  N 

angle  on  the  downstream  side  of  the  wave.  Now  the  change  in  deflection 
angle  between  D  and  C  is  given  by  the  relation  , 


r*D  -  ^  ^  + 


Thus  equations  (D2)  and  (D3 )  combine  to  yield 


dp  /db  dB  Ah’ 

PC‘pD  -  d&  nUc2  2+  &CxA  l 


But  the  compatibility  equations  (equations  (12)  and  (13))  combine  with 
equation  (Dl)  to  give  ^ 


H-4  =  -2®L  ( |l_  ^  zcx 

c  D  sin  2pN  \  dc2  '  dCi 


Equating  the  right-hand  members  of  this. and  the  previous  expression  and 
rearranging*  there  is  then  obtained 


dS/dCi 

5&/Sc2 


gy#N  _  d£ 

sin  2pw  dB  N 

’  2yNpN  |  dp 
sin  2pN  dB  N  , 


The  ratio  -  follows  from  the  sine  law;  however^  thus  the  ?  disturbance 
ACi  .  (  . 

strength  ratio  is  given  by  the  relation 


BS/dCi 

d5/6C2 


27npN  _  dp 
sin  2(3^  38  r 

27npN  dp 
in  2pN  d&  u 


,  * 

sin  (Pjj+bjj-a) 
sin  (I%-%hj) 


which  holds  for  both  ideal  and  calorically  imperfect  gas  flows 
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If  the  shock  wave  is  attached  to  the  airfoil  at  the  leading  edge 
(or  more  properly  if  MN>1),  this  expression  may  be  used  with  equa¬ 
tion  (23)  to  determine  the  surface  pressure  gradient  at  that  point.  A 
calculation  analogous  to  this  for  ideal  gas  flows  has  already  been  car¬ 
ried  out  by  Crocco,  Schaefer,  Munk  and  Prim  and  others  as  discussed  pre 
viously. 


NACA  TN  2646 


39 


REFERENCES 


1.  Ackeret,  J. :  Air  Forces  on  Airfoils  Moving  Faster  than  Sound. 

NACA  3M  317,  1925. 

2.  Busemann,  A.,  and  Walchner,  0.:  Airfoil  Characteristics  at  Super¬ 

sonic  Speeds.  British  R.T.P.  Trans.  1786.  (Forschung  auf  dem 
geheite  des  Ingenieurwesens,  vol.  4,  no.  2,  March/April  1933, 
pp.  87-92.) 

3.  Epstein,  Paul  S.:  On  the  Air  Resistance  of  Projectiles.  Proceed¬ 

ings  of  the  National  Academy  of  Sciences,  vol.  17,  1931, 
pp.  532-547. 

4.  Eggers,  A.  J.,  Jr.:  One -Dimensional  Flows  of  an  Imperfect  Diatomic 

Gas.  NACA  Rep.  959,  1950.  (Formerly  NACA  TN  l86l) 

5.  Ivey,  H.  Reese,  and  Cline,  Charles  W.:  Effect  of  Heat-Capacity  Bag 

on  the  Flow  Through  Oblique  Shock  Waves.  NACA  TN  2196,  1950. 

6.  Bethe,  H.  A.,  and  Teller,  E.:  Deviations  from  Thermal  Equilibrium 

in  Shock  Waves.  Ballistic  Research  Lab.  Rept.  No.  X-117. 

Aberdeen  Proving  Ground,  Aberdeen,  Maryland,  1945. 

7.  Crocco,  Luigi:  Singolarita  della  Corrente  Gassosa  Iperacustica 

Nell’ Intorno  di  una  Prora  a  Diedro.  L' Aerotecnica,  vol.  17, 
no.  6,  June  1937,  PP.  519-534. 

8.  Schaefer,  M. :  The  Relation  Between  Wall  Curvature  and  Shock  Front 

Curvature  in  Two-Dimensional  Gas  Flow.  USAF  Tech.  Rept. 
F-TS-1202-IA,  Jan.  1949. 

9.  Munk,  M.  M.,  and  Prim,  R.  C.:  Surface -Pres sure  Gradient  and  Shock- 

Front  Curvature  at  the  Edge  of  a  Plane  Ogive  with  Attached  Shock 
Front.  Jour.  Aero.  Sci.,  vol.  15,  no.  H,  Nov.  1948,  pp.  691-695. 

10.  Chapman,  Dean  R.:  Airfoil  Profiles  for  Minimum  Pressure  Drag  at 

Supersonic  Velocities  -  General  Analysis  With  Application  to 
Linearized  Supersonic  Flow.  NACA  TN  2264,  1951. 

11.  Isenberg,  J.  S.:  The  Method  of  Characteristics  in  Compressible 

Flow.  Part  I  (Steady  Supersonic  Flow)  USAF  Tech.  Rept. 
F-TR-1173A-ND,  Air  Materiel  Command,  Wright  Field,  Technical 
Intelligence  (Brown  University,  Graduate  Division  of  Applied 
Mathematics,  A-9-M  ll/l),  Dec.  1947* 


4o 


NACA  TN  2 646 


12.  Thomas ,  T.  Y.:  The  Determination  of  Pressure  on  Curved  Bodies 

behind  Shocks.  Communications  on  Pure  and  Applied  Mathematics, 
vol.  Ill,  no.  2,  June  1950,  pp.  103-132. 

13-  Noyes,  Robert  N. :  Prandtl-Meyer  Flow  for  a  Diatomic  Gas  of  Variable 
Specific  Heat.  NACA  TN  2125,  1950. 

14.  Tsien,  Hsue-shen:  Similarity  Laws  of  Hypersonic  Flows, 

Jour.  Math,  and  Phys.,  vol.  25,  no.  3,  Oct.  1946,  pp.  247-251. 

15*  Linnell,  Richard  D.:  Two-Dimensional  Airfoils  in  Hypersonic  Flows. 
Jour.  Aero.  Sci.,  vol.  16,  no.  1,  Jan.  1949,  pp.  22-30. 

16.  Lighthill,  M.  J. :  The  Flow  Behind  a  Stationary  Shock. 

Philosophical  Magazine,  ser.  7,  vol.  40,  Feb.  1949,  pp.  214-220. 

17*  Rand,  Robert  C«:  Prandtl-Meyer  Flow  Behind  a  Curved  Shock  Wave. 
Jour.  Math,  and  Phys.,  vol.  29,  no.  2,  July  1950,  pp.  124-132. 

18.  Busemann,  A.:  Flussigkeits-und  Garbewegung.  Handworterbuch  der 
Naturwissenschaften,  Zweite  Auflage,  Gustav  Fischer,  Jena,  1933, 
pp.  266-279- 

19-  Ivey,  H.  Reese,  Klunker,  E.  Bernard,  and  Bowen,  Edward  N. :  A  Method 
for  Determining  the  Aerodynamic  Characteristics  of  Two-  and  Three- 
Dimensional  Shapes  at  Hypersonic  Speeds.  NACA  TN  1613,  1948. 

20.  Meyer,  Th-:  The  Two-Dimensional  Phenomena  of  Motion  in  a  Gas  Flow¬ 
ing  at  Supersonic  Velocity.  Cornell  Research  Foundation,  Inc., 
Buffalo,  New  York,  Aug.  1945.  (Translation  from:  Mitteilungen 
uber  Forschungsarbeiten  auf  dem  Gebiete  des  Ingenieurvesens , 

Heft  62,  1908,  pp.  31-67.) 


TABLE  I.-  TABLE  OF  FUNCTIONS  FOR  SLENDER -AIRFOIL  METHOD 
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Figure  /.  -  Variation  with  deflection  angle  of  the  disturbance  strength  ratio  behind  an  oblique  shock 

wave  for  various  free -stream  Mach  numbers  (ys!.4). 


Pressure  coefficient,  C„  Pressure  coefficient. 
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(d)  Mo  =!0 

Figure  2.  -  continued. 


Pressure  coefficient,  Cp  Pressure  coefficient,  CP 
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Chordwise  station,  x/c 
(e)  Mo=/5 


Figure  2.  -  concluded. 


Chordwise  station,  x/c 
(f)  Mo  =  co 


Pressure  coefficient,  Cp  Pressure  coefficient. 
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Chordwise  station,  x/c 
(a)  Mo  =5 


Figure  3.  -  Pressure  distribution  on  iO  -percent-thick  biconvex  airfoil 

section  at  a=0? 


airfoil  section  at  a-  0° (rsi.4). 
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Figure  5.-  Range  of  applicability  of  shock -expansion  method  (y-1.4). 


Deflection  angle,  8,  deg  Deflection  angle,  8,  deg 
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Free- stream  Mach  number,  Mo 
(a)  Oblique  shock- wave  flows 


Free- stream  Mach  number,  Mo 
(b)  Expansion  flows 


Figure  6.  Accuracy  of  slender-airfoil  method  in  predicting  pressure 

coefficients  (y=l.4l 
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(c)  a  =  30.0  deg 

Figure  7.  ~  Pressure  distribution  on  10  -  percent- thick  biconvex 
airfoil  section  for  various  angles  of  attack  at  Mo-  to  (7s  14). 
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Angle  of  attack,  a,  deg 


(a)  Lift  coefficient. 


Angle  of  attack,  a,  deg 
(b)  Drag  coefficient. 


Angie  of  attack,  a,  deg 
(c)  Moment  coefficient. 


Figure  8. —  Variation  of  force  and  moment  coefficients  with  angle  of  oft 
for  !0-  percent  thick  biconvex  airfoil  section  at  M,  =  iO  (y=i.4). 
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Chordwise  station,  x/c 
(a)  a =  0° 


Chordwise  station,  x/c 


(b)  a  =  iO.O° 


Figure  II.  —  Pressure  distribution  on  lower  surface  of  10 -percent- thick  biconvex 
airfoil  section  for  various  angles  of  attack  at  Mo  =  10  and  To  -  500°  R. 


Pressure  coefficient,  C,  Pressure  coefficient,  Cp 
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Ca/oricai/y  imperfect  gas 
Shock-expansion  method 

- Slender- airfoil  method  Y=y 

Ideal  gas  * 

-Shock-expansion  method 

g 

Chordwise  station,  x/c 
(d)  a  -  30.0° 


Figure  //  -Concluded. 


fficient,  cm  Drag  coefficient,  cd  Lift  coefficient. 
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Angie  of  attack,  a,  deg 
(a)  Lift  coefficient. 
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Angie  of  attack,  a,  deg 
(b)  Drag  coefficient 


Caioricaiiy  imperfect  gas  / 

- Shock -expansion  method  _ ^ 

- Slender  airfoil  method  y='Xa  ^ 

Ideal  gas  _ 

- Shock - expansion ] 

method 
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Angie  of  attack,  a ,  deg 
(c)  Moment  coefficient 


Figure  12.- Variation  of  force  and  moment  coefficients  with  angle  of  attack 
for  iO~, percent- thick  biconvex  airfoil  section  at  M»  =  !0,  T,  =  50O°R. 
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Chordwise  station, 
(a)  Mo  =20 


(b)  Mo =30 


Figure  13.—  Pressure  distribution  on  iO-percent-thick  biconvex  airfoil  section 


at  a =0  and  1i=500°R 


'ariation  of  drag 


Pressure  coefficient. 
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Figure  15.  -  Pressure  distribution  on  iO- percent- thick  biconvex  airfoil  section 

at  Mo  =oo  and  a  =  0° 


me  18.  —  Accuracy  of  approximation  used  in  appendix  C  to  obtain  solution  for  Prandti  -  Meyer  flow  of  a  caioricaiiy 

imperfect  gas. 
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